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Introduction: We have developed a realistic, 

time-variable, one-dimensional, seasonal model for 
stratospheric photochemistry on Saturn using the Cal- 
tech/JPL KINETICS code [1,2,3]. The model ac- 
counts for variations in ultraviolet flux due to orbital 
position, solar-cycle variations, and ring- shadowing 
effects. The results for two Saturnian years, starting at 
L s = 0° in 1950 and running until the upcoming north- 
ern vernal equinox in 2009, are presented for numer- 
ous latitudes. The same two model years are run over 
and over again until the model convergences to make 
sure that high-altitude effects have had a chance to 
propagate down through the atmosphere. 

We use the SOLAR2000 model [4,5], in combina- 
tion with the spectra presented in [6], to predict the 
ultraviolet flux at any wavelength and any point in 
time during the simulation. Saturn’s orbital position 
during the simulation was taken from the ephemeris 
calculator at http://ssd.ipl.nasa.gov/horizons.html [7]. 
The photochemical model is derived from “Model C” 
of [8] and uses a hydrocarbon reaction list that has 
been extensively updated from that presented in [3]. 

Results: The effects of ring shadowing on the so- 
lar insolation are illustrated in Figure 1, where the 
derivation of the ring-shadowing equations is de- 
scribed in [9]. Several interesting effects can be noted 
from the figures. First, ring- shadowing effects are 
pronounced at low latitudes, become progressively less 
important at middle latitudes, and are unimportant at 
high latitudes. Second, daily insolation values are 
relatively constant with latitude in the summer hemi- 
sphere near summer solstice because the increase in 
the length of daylight hours with increasing latitude 
counteracts the increasing solar zenith angle. Latitude 
variations are more pronounced at other times of the 
year. Third, the fact that Saturn's orbit is eccentric and 
that perihelion occurs near southern summer solstice 
causes an asymmetry in the daily mean insolation be- 
tween the northern and southern hemispheres — 
higher daily mean insolation values are encountered 
during southern summer as opposed to northern sum- 
mer. By the same token, daily mean insolation values 
at low latitudes in the southern hemisphere exhibit 
more seasonal variation than low latitudes in the 
northern hemisphere. 

The latitudinal variation of the mixing ratios of 
C 2 H 6 and C 2 H 2 near southern summer solstice at dif- 
ferent pressure levels is shown in Figure 2. At the 
highest altitude levels (lowest pressures), the atmos- 
phere responds relatively quickly to insolation 


changes, and the hydrocarbon abundances have a lati- 
tude variation that resembles that of the solar insola- 
tion at L s = 273° (see also Fig. 1). The variation with 
latitude at lower altitudes (higher pressures) is more 
complicated; long diffusion time scales and phase lags 
affect the latitudinal profiles. Because vertical diffu- 
sion time scales are long in Saturn’s lower strato- 
sphere, the mixing-ratio profiles of photochemically 
stable hydrocarbons at pressures greater than ~0.1 
mbar in our model have a latitude variation that mim- 
ics the yearly average of the mean daily solar insola- 
tion, rather than mimicking the insolation variation 
from the current season. 
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Figure 1 . Daily mean solar insolation (W m' 2 ) incident on a 
unit horizontal surface at the top of Saturn’s atmosphere as a 
function of latitude and season ( L s ), assuming no ring shad- 
owing (top panel) or with the inclusion of ring shadowing 
(bottom panel). The shaded regions indicate polar night. L s 
= 0° at northern vernal equinox, 90° at northern summer 
solstice, 1 80° at northern autumnal equinox, etc. 

In Fig. 2, we also compare the models profiles with 
results from high-resolution, ground-based, infrared 
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observations using the TEXES grating spectrograph at 
the IRTF [10]. The shape of the observed C 2 H 2 distri- 
bution is similar to that of the model. However, the 
observed latitude distribution of C 2 H 6 does not com- 
pare well with the model. A quick examination of the 
relative photochemical lifetimes of the two species 
reveals a possible reason. The C 2 H 2 photochemical 
lifetime is shorter than the estimated vertical and me- 
ridional transport time scales and is similar to or 
shorter than a Saturnian season. The C 2 H 6 photo- 
chemical lifetime, on the other hand, is much greater 
than a Saturnian year and is longer than the estimated 
meridional and vertical transport time scales at the 
pressure levels at which the observations are sensitive. 
C 2 H 6 is therefore likely to be sensitive to transport 
effects and should act as a good tracer for atmospheric 
motions. 
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Figure 2. The latitudinal variation of the mixing ratios of 
C 2 H 2 (top) and C 2 H 6 (bottom) at L s = 273° (near southern 
summer solstice) at different pressure levels, as labeled. The 
model profiles (solid lines) are compared with results from 
ground-based infrared observations (diamonds) [10] that 
were acquired at similar L s . The open symbols in the model 
profiles represent models that have not yet converged. 

Conclusions: Realistic 1-D photochemical models 
that accurately track the time variability of incident 
solar ultraviolet radiation should be able to reliably 


predict (and explain) the latitudinal/seasonal behavior 
of relatively short-lived hydrocarbons like CH 3 , C 2 H 2 , 
C 2 H 4 , CH 3 C 2 H, and C 4 H 2 ; however, more observations 
are needed to test this claim. Longer-lived hydrocar- 
bons like C 2 H 6 and C 3 H 8 are more affected by atmos- 
pheric transport. Observations of the horizontal distri- 
bution of these long-lived species can help constrain 
stratospheric circulation patterns, and 2-D photo- 
chemical models will be needed to explain the ob- 
served latitudinal/seasonal behavior. Hydrocarbons 
located at pressures less than a few microbars will re- 
spond rapidly to changes in solar insolation; seasonal 
and latitudinal variations will be pronounced, and the 
ring shadow (and winter polar night) will strongly af- 
fect species abundances. In the middle and lower 
stratosphere, the latitudinal and seasonal distribution 
of hydrocarbons will be complicated and will depend 
on the relative magnitudes of the photochemical pro- 
duction and loss time scales, the vertical diffusion time 
scale, and the meridional transport time scale. Stable 
hydrocarbons residing at pressures greater than ~1 
mbar will be much less affected by latitude and season 
and will have latitudinal distributions that resemble 
that of the yearly average mean daily solar insolation. 

Cassini will map the hydrocarbon distributions on 
Saturn for the next several years, and these maps will 
provide a good test for our models, as well as good 
constraints for future 2-D models that include atmos- 
pheric transport. Ground-based infrared observations 
will continue to provide valuable information for sea- 
sons other than those encountered during the Cassini 
era and for altitudes that cannot be probed by the Cas- 
sini instruments. 
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